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electromagnetic  scattering  problems  involving  hailstones 
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hailstone  models  for  incident  wavelengths  of  1  and  10  cm. 

The  layered  hailstone  models  are  constructed  of  ice 

surrounded  by  a  layer  of  water  and  represent  melting 

hailstones.  Nonlayered  hailstone  models  represent 

nonmelting  hailstones  and  are  composed  entirely  of  ice. 

Hailstone  size,  shape,  and  structure  were  varied  to 

determine  their  effect  on  scattering  and  absorption 

characteristics.  The  results  indicate  that  resonance 

scattering  is  very  sensitive  to  hailstone  size,  shape,  and 

structure,  Rayleigh  scattering  is  sensitive  to  size  and 

structure  but  much  less  dependent  on  hailstone  shape. 
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SCATTERING  CALCULATIONS  FOR  VARIOUS  HAILSTONE  MODELS 
USING  THE  EXTENDED  BOUNDARY  CONDITION  METHOD 
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Overview 

This  thesis  determines  the  scattering  characteristics 
of  single,  resonance  sized,  axisymmetric  hailstone  models 
using  the  Extended  Boundary  Condition  Method  (EBCM).  This 

V 

theoretical  analysis  is  sponsored  by  the  Air  Force  Wright 
Aeronautical  Laboratories  (AFWAL)  at  Wright-Patterson  Air 
Force  Base.  This  chapter  presents  general  background 
material  and  a  brief  discussion  of  previous  work  on  the 
problem.  Information  concerning  the  approach  to  solving  the 
problem  is  also  presented. 

Background 

Secure  and  reliable  communications  systems  are  key 
elements  to  any  successful  Air  Force  operation;  success  or 
failure  of  a  mission  depends  on  accurate  communication  of 
information.  This  information  is  increasingly  being 
transmitted  via  millimeter  wave  satellite-to-earth 
communication  links.  Any  attenuation,  phase  shift,  or 
depolarization  of  the  transmitted  signals  prior  to  receive 
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must  be  considered  in  the  design  of  these  systems. 
Atmospheric  precipitation  in  the  form  of  rain  or  hail  can 
have  a  significant  effect  on  millimeter  wave  propagation. 
Scattering  by  rain  or  hail  also  plays  an  important  role  in 
radar  meteorology  (14:1). 

The  effect  of  precipitation  in  the  form  of  rain  has 
been  the  subject  of  extensive  research.  Considerably  less 
research  has  been  directed  toward  the  impact  of  hailstorms 
on  millimeter  wave  propagation  (14:1).  Scattering  effects 
of  hailstorms  can  be  determined  experimentally  or  predicted 
theoretically.  A  theoretical  investigation  of  hailstone 
scattering  involves  determining  the  individual 
electromagnetic  scattering  characteristics  of  the  hailstones 
and  modeling  the  distribution  of  the  hailstones  along  the 
propagation  path.  This  thesis  will  theoretically  examine 
the  individual  hailstone  scattering  problem.  In  the 
microwave  band  hailstones  are  often  resonance  sized  or  on 
the  order  of  a  wavelength  of  the  illuminating  signal.  A 
method  to  analyze  the  problem  must  be  selected  and  applied 
to  the  selected  hailstone  models. 

Scattering  characteristics  of  axisymmetric  objects  can 
be  theoretically  determined  using  several  techniques: 
integral  equation  formulations,  moment  methods,  and  finite 
element  methods  (6:227-323;  8:275-285;  9:202-214).  One  of 
the  most  accepted  integral  equation  techniques  is  the 
Extended  Boundary  Condition  Method  (EBCM)  introduced  by  P. 
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C.  Waterman  in  1965  (18:805-812).  The  EBCM  will  be  used  to 


handle  the  hailstone  models  considered  in  this  thesis. 

Problem 

This  thesis  will  use  the  EBCM  to  determine  the 
scattering  characteristics  of  single,  resonance  sized 
hailstone  models.  Wet  and  dry,  axisymmetric  models 
representing  melting  and  nonmelting  hailstones  will  be 
considered.  Sample  scattering  characteristics  of  the 
hailstone  models  will  be  presented. 

Summary  of  Current  Knowledge 

EBCM.  P.  C.  Waterman  proposed  the  EBCM  for  determining 
the  radar  cross  section  of  smooth,  perfectly  conducting 
bodies  illuminated  by  a  known  plane  wave.  The  EBCM  is  a  T- 
matrix  method  where  the  scattered  field  is  ultimately 
determined  through  the  use  of  a  transformation  matrix. 
Waterman  found  the  EBCM  to  be  well  suited  for  determining 
scattering  characteristics  of  smooth,  resonance  sized, 
axisymmetric  shapes  (18:811). 

Several  researchers  followed  up  on  Waterman's  original 
work  and  found  different  applications  for  the  method.  The 
EBCM  has  been  used  to  calculate  scattering  by  multilayered 
and  multiple,  dielectric  and  conducting  objects  including 
spheroids,  cones,  and  finite  cylinders  (3,5,11,15).  The 
next  chapter  contains  a  review  of  significant  contributions 
in  the  evolution  of  the  EBCM. 


Hail3tone  Models.  Hailstones  have  been  modeled  as 


several  different  shapes  to  study  their  scattering  affects. 
Among  these  are  spherical  models,  oblate  spheroidal  models, 
and  sphere-cone-spheroidal  models.  Layered  and  nonlayered 
models  representing  melting  and  nonmelting  hailstones  have 
been  considered.  The  nonlayered  models  are  composed 
entirely  of  ice  while  layered  hailstones  are  modeled  as  ice 
surrounded  by  a  film  or  layer  of  water  (1,4,15).  Pao  K. 

Wang  recently  proposed  a  three  parameter  hailstone 
representation  (16).  This  thesis  will  determine  scattering 
characteristics  of  Wang’s  hailstone  model.  The  next  chapter 
looks  at  the  various  hailstone  models  in  more  detail. 

Assumptions 

One  general  assumption  is  necessary  to  determine  the 
scattering  characteristics  of  the  hailstone  models.  For 
this  thesis  we  assume  the  hailstone  is  illuminated  by  a 
known  plane  wave.  This  assumption  is  necessary  so  the 
problem  can  be  analyzed  using  the  EBCM. 

Scope 

This  thesis  is  limited  to  axisymmetric ,  resonance  sized 
hailstone  models.  Nonlayered  hailstone  models  must  be 
homogeneous.  Although  layered  hailstone  models  will  be 
examined,  each  layer  must  itself  be  homogeneous.  The 
hailstone  models  considered  are  composed  of  pure  ice  and 
water  at  zero  degrees  Celsius. 


4 


Approach 

The  major  tasks  in  this  thesis  are  as  follows: 

(1)  Modify  an  existing  EBCM  program  to  determine 
layered  and  nonlayered  hailstone  model  scattering 
characteristics . 

(2)  Test  the  program  by  duplicating  previously 
published  results. 

(3)  Using  the  program,  determine  the  scattering 
characteristics  of  the  Wang  hailstone  models. 

Materials  and  Equipment 

An  EBCM  computer  program  to  be  modified  for  the 
hailstone  models  is  required  for  this  analysis.  Dr.  Peter 
Barber  at  the  Clarkson  College  of  Technology  graciously 
provided  his  EBCM  programs.  AFIT’s  CSC  and  ICC  computer 
systems  are  required  to  compile  and  execute  the  EBCM  code. 

Sequence  of  Presentation 

Chapter  II  will  review  the  literature  pertinent  to  this 
thesis.  Chapter  III  will  discuss  the  EBCM  theory  for 
layered  dielectric  scatterers  and  discuss  Wang’s  hailstone 
model.  Chapter  IV  will  present  the  scattering  results  for  a 
sample  of  hailstone  models.  Chapter  V  will  summarize  the 


effort  and  make  recommendations  for  future  research  efforts. 


II.  Review  of  Literature 


Introduction 

Electromagnetic  scattering  characteristics  of 
inhomogeneous,  axisymmetric  objects  can  be  theoretically 
determined  using  several  techniques:  integral  equation 
formulations,  moment  methods,  and  finite  element  methods 
(6:227-323;  8:275-285;  9:202-214).  Each  technique  has 
particular  advantages  and  limitations  with  respect  to 
scatter  size  and  shape.  One  of  the  most  accepted  integral 
equation  techniques  for  resonance  sized,  axisymmetric 
scatterers  is  the  Extended  Boundary  Condition  Method  (EBCM) 
introduced  by  P.  C.  Waterman  in  1965  (18:805-812).  This 
literature  review  examines  significant  contributions  in  the 
evolution  of  the  EBCM.  Various  hailstone  models  which  have 
been  studied  are  also  reviewed. 

Background 

Electromagnetic  scattering  problems  involve  an  incident 
wave  and  a  scattering  object.  The  incident  wave  impinges 
upon  the  object  and  sets  up  currents  in  and  on  the  object. 
These  induced  currents  create  an  electromagnetic  field  which 
radiates  in  all  directions.  The  nature  of  the  incident  wave 
and  the  geometry  and  composition  of  the  scattering  object 
determine  the  characteristics  of  the  scattered 
electromagnetic  field.  Maxwell’s  equations  and  appropriate 


boundary  conditions  can  be  used  to  relate  the  induced 
current  and  scattered  field  to  the  incident  wave  and 
scattering  body  (7:52). 

Scattering  Regimes.  Electromagnetic  wave  scattering 
problems  are  characterized  by  the  relationship  between 
incident  wavelength  and  scatterer  size.  The  spectrum  of 
scattering  problems  is  divided  into  three  regimes:  the 
Rayleigh  region  where  incident  wavelength  is  much  larger 
than  the  body  size,  the  resonance  region  where  incident 
wavelength  is  on  the  order  of  body  size  and  the  optics 
region  where  incident  wavelength  is  much  smaller  than  the 
body  size  (7:53).  In  general,  an  analytic  technique  for 
solving  scattering  problems  is  applicable  to  one  of  the 
three  regimes.  The  EBCM  is  most  applicable  to  scattering 
problems  in  the  resonance  region  (3:2864). 

Resonance  Region  Scattering.  Figure  1  shows  a  typical 
scattering  problem.  The  scatterer,  S,  is  illuminated  by  a 
plane  wave,  E^  and  ,  traveling  in  direction  k.  Es  and  Hs 
are  the  resulting  scattered  fields.  For  a  resonance  sized 
scatterer,  the  phase  of  the  incident  field  changes  over  the 
length  of  the  scatterer.  Analyzing  a  resonance  region 
scattering  problem  requires  considering  the  affect  of  each 
part  of  the  scatterer  to  every  other  part  (7:57).  Rayleigh 
and  optics  approximations  can  not  be  used  since  the  complete 
wave  nature  of  the  incident  field  must  be  considered  in  the 
scattering  problem  (3:2864). 
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Figure  1.  Electromagnetic  Scattering 

Origin  of  the  EBCM 

P.  C.  Waterman  proposed  a  new  method,  the  EBCM,  for 
determining  the  radar  cross  section  and  other  scattering 
quantities  for  smooth,  perfectly  conducting  bodies 
illuminated  by  an  incident  plane  wave  (18:805).  In  the 
EBCM,  known  and  unknown  fields  are  expanded  into  infinite 
series  of  spherical  vector  wave  functions.  Numerically,  the 
infinite  series  expansions  are  truncated  at  some  value 
determined  by  convergence  testing.  The  unknown  expansion 
coefficients  are  solved  for  in  terms  of  the  known 
coefficients  via  a  transformation  or  T-matrix.  Waterman 
implemented  the  method  on  a  computer  and  determined  the 
radar  cross  section  (RCS)  of  a  sphere-cone-sphere.  He  found 
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the  EBCM  to  be  well  suited  for  handling  smooth  axi symmetric 
shapes.  Computation  of  surface  current,  scattered  field, 
and  differential  scattering  cross  section  were  also 
accomplished  with  little  additional  computer  time  (18:811). 
Waterman  noted  no  deterioration  of  results  as  frequency 
increased  with  respect  to  scatterer  size  and  that  no 
numerical  difficulties  were  observed  for  resonance  sized 
objects  (18:811).  Although  resonance  was  not  a  problem, 
Waterman  observed  "poor  numerical  convergence  of  the 
truncation  procedure"  when  the  method  was  applied  to 
elongated  shapes  with  a  length-to-width  greater  than  two 
( 18:811 ) . 

In  a  later  article,  Waterman  expanded  his  method  to 
compute  the  RCS  of  two  types  of  conducting  objects  with 
edges:  cone-spheres  and  right  cylinders.  The  theoretical 
results  from  the  cone-sphere  and  cylinder  compared  well  with 
experimental  measurements.  Perfectly  conducting,  prolate 
and  oblate  spheroids  were  also  examined  and  compared  to 
geometrical-optics  solutions.  Agreement  of  the  two  methods 
was  not  good  for  resonance  sized  spheroids.  This  was 
expected  because  geometrical-optics  methods  are  generally 
not  applicable  to  problems  in  the  resonance  region.  The  two 
methods  compared  well  for  larger  spheroids  (19:825-839). 

Multilayered  and  Dielectric  Scatterers 


Waterman’s  work  had  been  limited  to  homogeneous 
perfectly  conducting  bodies.  Peterson  and  Strom  showed 


Waterman’s  method  could  be  extended  to  apply  to  multilayered 
and  multiple,  dielectric  scatterers.  Peterson  and  Strom’s 
work  was  mostly  theoretical  but  some  numerical  computation 
was  accomplished.  Resonance  sized  dielectric  spheres  were 
examined.  Also,  the  RCS  of  a  perfectly  conducting  sphere 
embedded  in  a  concentric  dielectric  sphere  was  computed. 

The  EBCM  produced  exact  solutions  for  scattering  by  spheres 
(  11:2670-2684)  . 

Bringi  and  Seliga  followed  Peterson  and  Strom’s  work  on 
multilayered,  axisymmetric  scatterers  (5:575-580).  Unlike 
Peterson  and  Strom  who  used  the  Poincare-Huygens  principle, 
Bringi  and  Seliga  derived  the  EBCM  equations  for 
multilayered  scatterers  using  the  equivalence  principle. 
Bringi  and  Seliga  used  the  method  to  compute  scattering 
characteristics  of  dielectric  spheres  or  conductors  embedded 
in  a  spherical  dielectric  shell.  The  results  compared  well 
to  extended  Mie  theory  calculations.  Bringi  and  Seliga 
thought  the  method  would  be  well  suited  to  examine  the 
problem  of  scattering  by  melting  hailstones  (5:580). 

Independent  of  Peterson  and  Strom,  Barber  and  Yeh 
solved  axisymmetric,  dielectric  scattering  problems  using 
the  EBCM  (3:2864-2872).  The  scattering  objects  considered 
were  resonance  sized,  dielectric  spheres,  prolate  and  oblate 
spheroids,  and  finite  cylinders  (3:2864-2872).  In  a  later 
work,  Barber  proposed  a  method  for  extending  the  region  of 
applicability  of  the  EBCM  for  highly  elongated  or  high  loss 
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objects.  Barber’s  proposal  consisted  of  making  a  set  of 
calculations  for  a  series  of  reduced  dielectric  constants 


where  good  results  were  found.  The  desired  characteristics 
could  then  be  extrapolated  from  a  resulting  curve  for  the 
dielectric  constant  of  interest  (2:380). 

Other  independent  work  on  Waterman's  method  was 
conducted  by  Warner  and  Hizal  (17:921-930).  Warner  and 
Hizal  were  interested  in  applying  the  EBCM  to  determine 
attenuation  and  depolarization  effects  of  rain  on  microwave 
transmissions  (17:921).  Single  raindrops  were  modeled  as 
homogeneous,  dielectric  oblate  spheroids  and  results  from 
the  EBCM  compared  well  with  previously  published  work. 

Warner  and  Hizal  concluded  that  Waterman’s  method  could  be 
applied  to  hailstone  models  as  well  (17:928). 

Hailstone  Models 

Hailstones  of  various  size  have  been  modeled  as 
axisymmetric ,  geometrical  shapes.  Among  these  are  spherical 
models,  oblate  spheroidal  models,  and  sphere-cone-spheroidal 
models.  Pao  K.  Wang  recently  proposed  a  three  parameter 
hailstone  representation.  His  model  is  different  than 
previous  models  in  that  size  and  shape  distributions  can  be 
described  more  quantitatively  (16:1062).  Figure  2  shows  the 
general  shapes  of  different  hailstone  models. 

Spherical  Hailstone  Model.  Battan,  Browning,  and 
Herman  modeled  hailstones  as  spheres  and  determined 
scattering  characteristics  using  the  extended  Mie  theory. 


\m*w$ 


Wet  and  dry  hailstone  models  were  considered.  The  wet 
model,  representing  a  melting  hailstone,  was  a  sphere  of  ice 
surrounded  by  a  slightly  larger  diameter  shell  of  water  (4). 
Battan  Browning,  and  Herman  were  the  first  to  examine  the 
hailstone  scattering  problem  and  their  results  provided  a 
basis  for  comparison  of  future  work. 


Figure  2.  Hailstone  Models 


theory  used  by  Battan,  Browning,  and  Herman  can  be  used  to 
solve  problems  involving  concentric  spheres.  Other  methods 
must  be  used  for  aspherical  shapes.  Wang  and  Barber  used 
the  EBCM  to  determine  scattering  characteristics  of  oblate 
spheroidal  models  of  hailstones.  They  compared  layered 
oblate  spheroidal  models  of  varying  size  and  axial  ratio 
with  corresponding  equivolume  spheres.  Although  absorption 
behavior  was  similar,  the  oblate  spheroidal  models  generally 
scattered  more  incident  radiation  than  an  equivolume  sphere. 
For  axial  incidence,  backscattering  increased  significantly 
with  increased  axial  ratio.  Wang  and  Barber  concluded  that 
careful  modeling  of  layered  physical  objects  is  important  to 
obtain  accurate  theoretical  calculations  (15:1190-1197). 

Sphere-Cone-Spheroidal  Hailstone  Model.  Aydin  and 
Seliga  modeled  hailstones  as  sphere-cone-spheres  and  sphere- 
cone-oblate  spheroids.  They  used  the  EBCM  for  their 
backscatter  analysis  and  also  examined  the  oblate  spheroid 
hailstone  model  used  by  Wang  and  Barber.  Layered  and 
homogeneous  models  representing  melting  and  nonmelting 
hailstones  were  considered  at  3  GHz.  Effects  due  to  various 
water  coating  thicknesses  were  examined  (1:58-66).  The 
scattering  characteristics  of  the  sphere-cone-spheroidal 
hailstone  models  were  compared  to  equivolume  spherical 
models.  Backscattering  characteristics  varied  with 


hailstone  size,  shape  and  composition. 


Wang's  Hailstone  Representation.  Pao  K.  Wang’s  three 
parameter  hailstone  representation  was  used  to  analyze  the 
size  and  shape  distribution  of  a  sample  of  679  hailstones. 
The  hailstones  were  collected  during  a  hailstorm  in 
Colorado.  The  analysis  revealed  simple  relations  between 
the  number  concentrations  of  hailstones  and  each  of  the 
three  hailstone  parameters.  This  study  was  concerned  with 
size  and  shape  characteristics  of  hailstones  (16:1062-1070). 
The  Wang  hailstone  model  will  be  discussed  in  more  detail  in 
the  next  chapter.  Theoretical  scattering  characteristics  of 
Wang’s  hailstone  model  have  not  been  determined. 

Summary 

Several  applications  and  their  contribution  to  the 
evolution  of  the  EBCM  have  been  reviewed.  The  method  has 
been  used  to  determine  scattering  and  absorption 
characteristics  of  resonance  sized,  perfectly  conducting  and 
dielectric,  axisymmetric  bodies.  Multilayered  bodies  and 
objects  larger  than  incident  wavelength  have  been  examined 
using  the  EBCM.  Scattering  characteristics  of  spherical  and 
oblate  spheroidal  models  of  wet  and  dry  hailstones  have  been 
determined.  Scattering  characterization  of  Wang’s  hailstone 
representation  have  not  been  reported. 

This  chapter  presented  a  review  of  literature 
pertaining  to  the  EBCM  and  hailstone  models.  The  next 
chapter  presents  the  theory  of  the  EBCM  for  layered, 


t 1  <  i1 


III.  Theory 

Overview 

This  chapter  summarizes  the  theory  of  the  EBCM  for  a 
two-layered,  axisymmetric ,  dielectric  scatterer .  The 
notation  and  most  equations  in  this  summary  are  taken  from 
the  EBCM  theory  presented  by  Wang  and  Barber  (15).  Consult 
Peterson  and  Strom  (11)  for  a  detailed  discussion  of  the 
method  for  multilayered  scatterers.  This  chapter  also 
discusses  the  size  and  shape  characteristics  of  hailstones 
using  Wang’s  three  parameter  hailstone  representation. 

The  Scattering  Problem 

Figure  3  is  a  model  of  the  problem  of  electromagnetic 
scattering  by  a  two-layered  object.  For  this  problem  there 


E3  =  E1  +  E3,  H3  =  H1  +  H3 


Figure  3.  Layered  Scattering  Problem 
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are  three  regions  separated  by  object  surfaces  and  S2* 

E*,  H*  ;  E^,  H^;  and  E^,  are  the  fields  in  regions  1,  2 
and  3  respectively.  E*,  H*-  is  the  known  incident  field  and 
Es ,  Ha  is  the  unknown  scattered  field.  As  shown  in  the 
figure,  the  total  field  in  region  3,  (E^,  §3 ) ,  is  the  sum  of 
the  incident  and  scattered  fields.  Each  of  the  three 
regions  of  the  problem  has  the  permeability  of  free  space, 

M0 ,  and  a  complex  permittivity,  €^,  where  i  =  1,  2,  or  3. 
Each  region  also  has  an  associated  complex  wave  number  kj. 

ki  =  (2n/X0) (€i/€0)1/2  (1) 

where  Xq  is  the  free-space  wavelength.  Outward  directed 
normal  unit  vectors  to  surfaces  Sj  and  S£  are  shown  as  rtj 
and  2 • 

Our  goal  is  to  determine  the  scattering  characteristics 
of  the  object.  To  achieve  this  goal  we  will  solve  for  the 
scattered  field  in  terras  of  the  known  incident  field  and 
scatterer  properties.  The  fields  will  be  expanded  using 
spherical  vector  wave  functions.  The  expansion  coefficients 
of  the  incident  field  are  known  and  all  other  field 
expansion  coefficients  are  unknown.  The  solution  ultimately 
relates  the  scattered  field  expansion  coefficients  to  the 
known  incident  field  expansion  coefficients  via  a 
transformation  or  T-matrix.  To  find  our  solution  using  the 


EBCM,  we  will  make  the  following  assumptions:  the  object  is 
axisymmetric  and  the  complex  permittivity  within  each  region 
is  constant  throughout  that  region.  We  also  assume  an 
exp(-jwt)  time  variation. 

Theoretical  Development 

The  equivalence  theorem  permits  us  to  reduce  the 

original  problem  of  Figure  3  to  the  three  subproblems  shown 

in  Figure  4.  For  each  subproblem,  the  surface  currents  J 

and  M  produce  the  actual  field  in  the  region  of  interest  and 

a  null  field  elsewhere.  The  plus  (+)  and  minus  (-) 

subscripts  indicate  external  and  internal  field  quantities, 

respectively,  evaluated  at  the  surface. 

First,  we  examine  the  subproblem  shown  in  Figure  4a. 

The  scattered  field  external  to  S2  is  due  to  surface 
-*•2  -*-2 

currents  J+  and  M+  and  can  be  found  using  vector  potentials 
A  and  F . 


Es  =  -  V*  F  -  — —  (  V  x  V  x  A) 
jw€0 


(2a) 


Hs  =  V  *  A  -  _  (V*  V  *  F ) 
JWUQ 


(  2b) 


where 


As  _ 


1 


f  J+e*p( J^3 \t  ~  f  1 ) 


H  Q  * 


(  \ 
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S  -  1  f  M+exp( jk3 j*  -  *’|) 

-  k  js-2 — 


(3b) 


&nd 


jJ  =  ft2  x  hJ  (4a) 

M+  =  -  !t2  *  E+  (4b) 

The  position  vectors  *  and  *’  in  Equation  2  indicate 
observation  and  source  points  respectively. 

The  total  field  external  to  S£  in  Figure  2a  is  the  sum 
of  the  incident  and  scattered  fields.  So,  for  *  outside  S2 
we  have 

E 3(*)  =  E i(t)  +  E 3(t)  (5a) 

A  null  field  exists  within  S2,  so  for  *  inside  S2  we  have 

0  =  EM*)  +  E3  ( * )  (5b) 

Substituting  Equations  2  through  4  into  Equation  5  yields 
the  following  set  of  integral  equations 
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EM*) 


+  Vx  j 


S2 


(*2  x  E+ ) g ( k3R )ds ' 


-  Vx  Vx  /  -±-  ( *2  x  Hj)g(k3R)ds’ 

J S2Jw€3 


J 


•I 


E^(*)  ;  *  outside  S2 
0  ;  *  inside  S2 


( 6e ) 
(6b) 


where  g(k^R)  is  the  free-space  Green’s  function  for  region  i 


g(kiR)  = 


exp( jkjR) 
4kR 


(7a) 


R  =  I*  -  *’l 


(7b) 


Now  we  look  at  the  subproblem  shown  in  Figure  4b.  In 
this  subproblem,  the  surface  currents  produce  a  field  in  the 
region  between  surfaces  and  S2  and  a  null  field  outside 
S2  and  inside  Sj.  The  following  set  of  integral  equations 
results  from  analysis  of  this  subproblem 


x  E?)g(k2R)ds’ 
x  E+)g(k2R)ds 


Vx  Vx  /  -L-  (-*2  x  H?  )  g  ( k2R )  ds  ' 

J S2jw€2 


Vx  Vx  /  -1-  (*!  X  )  g  (  k2R )  ds 

jw€2 


E  (*)  ;  *  between  Si  and  S2  (8a) 

0  ;  t  outside  S2  and  inside  Sj  (8b) 


We  now  have  integral  equation  representations  of  the 
subproblems  of  Figure  4a  and  4b.  Similar  equations  can  be 


written  for  the  subproblem  of  Figure  4c  but  are  not  necessary  in 
this  analysis.  Boundary  conditions  exist  at  object  surfaces 
and  S£  enabling  us  to  relate  Equations  7  and  8  and  ultimately 
solve  the  scattering  problem.  At  object  surfaces  Sj  and  S£  the 
tangential  fields  must  be  continuous.  Thus,  we  have  the 
following  boundary  conditions: 

*  E+  =  *  E-  (9a) 

1^2  x  H+  =  x  H-  (9b) 

ltl  x  E+  =  x  e!  ( 9c ) 

ft}  x  H+  :  ft}  x  5!  ( 9d ) 


Substituting  Equations  9a  and  9b  into  Equation  6  and 
Equations  9b  and  9c  into  Equation  8  we  have 


Ei(f)  +  Vx  I  <!t2  x  E?) «S(k3R)ds 

J  S2 


-  V*  Vx  /  -±-  <rt2  x  H?)  •G(k3R)ds’ 

J S2jw€3 


•I 


E  (t)  ;  t  outside  S2 
0  ;  t  inside  S2 


( 10a) 
(  10b) 


C 


c 
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-rtf* 


Vx 


/  ( -ft2  x  E?)  •ff(k2R)ds'  -  Vx  Vx  /  -L_  (-ft2  x  H?).5(k2R)ds 

J s2  y s2jw€2 


Vx  /  <  X  e!)  •ff{k2R)ds'  -  Vx  Vx  /  -A-  (*!  X  h!  )  •2f(k2R 

ysj  Js1J"€2 


)ds 


■{ 


E  *£)  ;  1*  between  Sj  and  S2  (11a) 

0  ;  £  outside  S2  and  inside  Sj  (lib) 


The  free-space  Green's  dyadic,  2f(kiR),  was  used  in 
Equations  10  and  11.  This  is  necessary  to  numerically  solve  the 
integral  equations  for  the  scattered  field.  The  dyadic  Green’s 
function  and  the  fields  in  Equations  10  and  11  will  be  expanded 
in  spherical  vector  wave  functions  M  and  N  (13:415). 


>  mL(W)  =  Vx  t  C08;m*!  Pn(cose)2P(kr)  (12a) 

’  sin(m«) 

Nomn<*>  =  r  Vx  Mg^t)  (12b) 

k 

► 

where  pJJ(cos®)  is  an  associated  Legendre  function,  a  is  even 
or  odd,  and  zg(kr)  is  an  appropriate  spherical  Bessel 
}  function. 


5P(kr)  =  f  Jn<^) 

\  jn<kr)  +  Jnn( 


;  p  =  1 

kr )  ;  p  =  3 


(13) 


The  dyadic  Green’s  function  expansion  is  given  as 


G(  kR )  =  —  T'Dv[Mj(k£>)Mi(k*<)  ♦  Ny(  k£>  )  Ny(  k?<  )  ]  (14) 

-*r 


where  *>  is  the  greater  of  f  or  t’ ,  *<  is  the  lesser  of 
t  or  t’ ,  and  Dv  is  a  normalization  constant  (10:1781).  The 
summation  over  v  is  a  combined  summation: 


E  =  23  E  E 


a rg  n= 1  m=0 


= 


_  €m( 2n  +  1 ) (n  -  m) ! 


4n(n  +  l)(n  +  m)! 


•m 


■I 


m  =  0 
m  >  0 


The  fields  in  Equations  10  and  11  are  also  expanded  in 
of  M  and  N. 


EM*)  =  EDv[avMi(k3*)  +  bvMi(k3*)J 

v 

E*(t)  =  E 4Dv[fvMj(k3*)  +  gvNv( k3*) ] 

v 

E2<*)  =  ElMTvMi(k2*)  +  6uNi(k2*) 

v 

+  avM2(k2*)  +  flvNj(k2*)l 
E*(*)  =  E  DvfCvM^ki*)  +  d VN y ( k x * )  ] 

V 

H2(t)  =  -j(€2/u0)1/2  E  Dv[TvNi,(k2*)  +  6vMl(k2*) 

V 

avN2(k2t)  +  8vM„(k 2*)1 


(15a) 

(  15b) 

(  15c) 

terms 

(16a) 

(  16b) 

(16c) 

(  1 6d ) 


(  1 6e  ) 


H^t)  =  -j(€i/M0)1/2  EDv[cvNl(ki?)  +  dvM^kif)]  (16f) 


The  following  relations  were  used  to  determine  Equations  16e 
and  16f  from  Equations  16c  and  16d  respectively  (13:715): 


r  V  x 


(17a) 


i  Vx  N 
k 


(  17b) 


-  Vx  M 
k 


(  17c) 


Substituting  Equations  16  into  Equations  10  and  11  and 
truncating  the  summations  over  v  at  some  upper  index  N 
enables  us  to  numerically  solve  for  the  scattered  field 
expansion  coefficients,  f  and  g,  in  terms  of  incident  field 
coefficients,  a  and  b. 


f  a/4 

g  =  - CT3  b/4 


[T]  is  the  2N  x  2N  T-matrix  containing  the  scatterer 
information.  It  is  given  by 


(T]  =  { CQ2 1 1  ~ 


*  { [Q2  1 3  -  [Q23HD][q}1][Qi1]-1(D]'1}'1  (19) 
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where  the  following  notation  has  been  used: 


pq 


pq 

Hi  3Wv 


pq 

[Ji  lWv 


pq 

[Ki  ]uv 


pq 

tLi  l*iv 


Ur.) 


1/2 


t  pq 

♦  (Cr,)1^ 

«.r*  * 

Ki 

(  20a) 

_ir 

*  (€ri)l/2L?q 

^  ♦  Uri|l/2K?\ 

.2 

ki  +  l  j 

f  Iti-  M^k i  +  1t) 

x  My(ki^)ds 

(20b) 

n  J 

Si 

.2 

ki  +  l  i 

f  J V  Mj(ki  +  1f) 

x  Ny(k^^)ds 

(20c) 

*  J 

fSi 

.2 

ki  +  l  1 

f  rti-  NjS(ki+1t) 

x  My(ki*)ds 

(  20d ) 

*  J 

Si 

ki  +  lj 

f  IV  N  ^  ( k  i  + 1  ^ ) 

x  Nyjki'fjds 

(20e) 

*  J 

fSi 

Ui)1/2/<n*:>1/2 

( 20f ) 

Normalization  constants  are  contained  within  the  diagonal 
matrix  [D] . 

Different  scattering  quantities  can  be  calculated  after 
determining  scattered  field  coefficients  f  and  g.  The 
vector-far-f ield  scattered  field  amplitude  F  can  be  found 
from  the  following  relationship: 


E3(k t)  =  F(es>03;ei,0i)ex?(-J-k_rj  ;  as  kr  - 


(21) 


where  (ea,*a)  is  the  scattered  field  direction  and  ( ®i » *i ) 


is  the  incident  field  direction.  Other  calculable 


quantities  include  differential  scattering  cross  section 
(DSCS),  extinction  cross  section  ( °ext ) >  scattering  cross 
section  ( osc )  and  absorption  cross  section  ( oabs ) •  These 
quantities  are  obtainable  from 


DSCS  = 


(22a) 


°ext  =  <  ®i  >  *i  J  ®i  «  *i  )  1 

k 


(22b) 


'sc 


16* 

~k2 


£Dv[  fy  2  +  *v  21 


(22c) 


°abs  =  °ext  -  °sc 


( 22d ) 


where  e  is  the  polarization  vector  of  the  incident  field. 
Other  quantities,  such  as  radar  cross  section,  can  also  be 
calculated  from  scattered  field  coefficients  f  and  g. 

Wang’s  Hailstone  Model 

The  next  chapter  will  discuss  the  scattering 
characteristics  of  Wang’s  three  parameter  hailstone 
representation.  This  section  provides  some  physical  insight 
for  the  model  itself. 

The  following  mathematical  expression  describes  a  curve 
that  becomes  the  surface  of  Wang’s  hailstone  model  when 
rotated  about  the  z-axis  (16): 


2 


From  Equation  23  and  Figure  5  we  can  see  how  parameters  a, 
c,  and  x  affect  the  shape  of  the  hailstone.  Increasing  a 
while  holding  c  and  x  constant  causes  the  hailstone  to 
become  broader  along  the  x-axis.  Increasing  c  causes 
elongation  in  the  z  direction.  Parameter  x  is  more 
difficult  to  understand.  Figure  5  was  drawn  with  shape 
factor  x  =  1.  As  x  increases,  the  hailstone  rapidly  becomes 
more  spheroidal.  Figures  6  and  7  show  hailstones  with  x=2 
and  x=10  respectively.  Parameters  a  and  c  are  constant  for 
Figures  5  through  7 . 


Figure  6.  Hailstone  Model  x=2 


Figure  7.  Hailstone  Model  x=10 

A  sample  consisting  of  679  hailstones  was  collected  and 
analyzed  using  Wang’s  hailstone  representation  (16). 

Figures  8,  9,  and  10  show  the  distribution  of  parameters  a, 
c  and  X,  respectively,  for  the  sample.  The  (a)  parameter 
varied  in  size  from  about  0.15  to  0.65  cm  for  the  sample. 
Approximately  one-third  of  the  hailstones  had  an  (a) 
parameter  between  0.25  -  0.3  cm.  Parameter  (c)  varied  from 
0.3  to  0.9  cm  with  the  interval  0.5  to  0.6  cm  containing 
greater  than  one-third  of  the  sample.  Shape  parameter,  (X), 
varied  from  1  to  10.  Half  the  hailstone  sample  had  a  shape 
parameter  (x)  between  1  and  2  (16). 
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Figure  10.  The  X-distribution  (16:1066) 


The  melting  hailstones  are  modeled  as  ice  layered  with 
varying  thicknesses  of  water.  The  outer  dimensions  of  the 
layered  models  are  held  constant  as  thickness  (t)  is  varied. 
Water  layer  thickness  (t)  is  shown  in  Figure  11.  No 
experimental  data  are  available  for  parameter  (t). 

Scattering  by  melting  hailstones  will  be  considered  in  the 
next  chapter. 


Figure  11.  Layered  Hailstone  Model 


This  chapter  summarized  the  EBCM  theory  for  a  layered 
scatterer.  The  T-matrix  for  determining  scattered  field 
coefficients  from  incident  field  coefficients  was  presented 
and  calculable  scattering  quantities  were  discussed.  Size 
and  shape  characteristics  of  Wang's  hailstone  model  were 
also  presented.  The  next  chapter  examines  how  scattering 
characteristics  of  the  hailstone  models  change  as  the  model 
parameters  are  varied. 


IV.  Results 


Overview 

This  chapter  presents  and  discusses  scattering 
characteristics  of  Wang's  hailstone  model.  An  existing  EBCM 
computer  program  for  determining  scattering  by  spheroids  was 
modified  to  obtain  the  scattering  plots  for  Wang’s  model. 
Scattering  information  is  presented  for  a  sample  of 
hailstone  models  at  3  and  30  GHz. 

Program  Modification 

Dr.  Peter  Barber  from  the  Department  of  Electrical  and 
Computer  Engineering  at  the  Clarkson  College  of  Technology 
provided  computer  code  for  calculating  scattering  by  layered 
and  nonlayered  dielectric  spheroids  using  the  EBCM. 
Documentation  and  sample  computer  runs  were  provided  with 
the  FORTRAN  code.  The  sample  computer  runs  included  input 
and  output  for  a  nonlayered  prolate  spheroid. 

The  program  was  compiled  on  the  AFIT  Hercules  computer 
system.  The  Hercules  is  a  Digital  Equipment  Corporation 
VAX-11/785  supercomputer  running  the  VMS  operating  system. 
The  prolate  spheroid  test  cases  were  run  and  verified. 

The  code  provided  by  Dr.  Barber  determines  scattering 
by  spheroidal  dielectric  objects.  Spheroids  are 
axisymmetric  as  well  as  mirror  symmetric.  Top  and  bottom 
halves  of  mirror  symmetric  objects  have  the  same  physical 
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dimensions.  Figure  12  illustrates  mirror  and  nonmirror 
symmetry.  The  Wang  hailstone  model  is  axisymmetric  but 
generally  not  mirror  symmetric.  Dr.  Barber’s  programs  were 
modified  for  nonmirror  symmetric  objects  to  determine  the 
scattering  characteristics  of  the  Wang  model. 


Figure  12.  Mirror  and  Nonmirror  Symmetry 


Several  test  cases  were  run  to  verify  the  nonmirror 
symmetric  code.  The  original  prolate  spheroid  test  case 
produced  the  same  result  and  Aydin  and  Seliga’s  results  were 
duplicated.  Aydin  and  Seliga  examined  differential  radar 
scattering  properties  of  several  hailstone  models  (1:58-66). 
The  hailstone  models  examined  include  layered  and  nonlayered 
oblate  spheroids,  nonlayered  sphere-cone-spheres,  and 
nonlayered  sphere-cone-oblate  spheroids.  The  sphere-cone- 


sphere  and  sphere-cone-oblate  spheroidal  hailstone  models 
are  nonmirror  symmetric.  A  sample  of  Aydin  and  Seliga’s 
work  was  duplicated  and  is  presented  in  Appendix  A. 

The  programs  provided  by  Dr.  Barber  were  capable  of 
computing  up  to  20-by-20  T-matrices.  This  corresponds  to 
truncating  the  field  expansions  at  a  maximum  N  value  of  10. 
Larger  N  values  were  needed  to  compute  the  T-matrices  for 
the  Wang  hailstone  models.  In  general,  as  scatterer  size, 
asymmetry  and  dielectric  constant  increase  the  required  N 
and  resulting  T-matrix  size  increase.  The  programs  were 
modified  to  compute  T-matrices  for  a  maximum  N  value  of  40. 
Complex  double  precision  was  also  necessary  to  compute  the 
spherical  Bessel  functions.  Computing  time  increases 
significantly  with  increasing  N.  The  AFIT  ELXSI  6400  was 
used  for  N  values  greater  than  20  due  to  its  faster  and  more 
precise  fortran  code  execution.  A  brief  description  of  the 
final  code  used  in  this  thesis  is  included  in  Appendix  B. 

Coordinate  System 

Figure  13  shows  the  coordinate  system  used  in  this 
analysis.  The  incident  wave  is  linearly  polarized  and 
propagates  in  the  +z-direction .  In  general,  the 
polarization  vector  of  the  incident  wave  lies  in  the 
xy-plane  and  is  specified  by  magnitude,  A,  and  angle,  v>, 
from  the  x-axis.  For  this  analysis,  the  polarization  vector 
of  the  incident  wave  is  fixed  with  A  =  1  and  v>  ~  0  degrees. 


The  orientation  of  the  scatterer  is  specified  by  angles 


©  and  «  where  S  is  the  axis  of  symmetry  of  the  scatterer. 

The  top  of  the  scatterer  is  in  the  +S-direction .  The  center 
of  the  scatterer  is  located  at  the  origin  of  the  coordinate 
system.  The  scatterer  is  fixed  at  ©  =  «  =  90  degrees  for 
horizontal  polarization  calculations  and  ©  =  90  degrees, 

*  s  0  degrees  for  vertical  polarization  calculations. 


Figure  13.  Coordinate  System 


The  scattered  field  is  evaluated  in  xz-plane  with  0 
degrees  in  the  forward  or  +z-direction  and  180  degrees  in 
the  backscattered  or  -z-direction .  The  scattering  angle  is 
equal  to  90  degrees  at  the  +x-axis. 


Graphical  Results 


Scattering  characteristics  of  various  hailstone  models 
were  calculated  at  two  frequencies:  3  GHz  and  30  GHz.  The 
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results  are  presented  in  graphical  form  in  this  section. 

The  hailstone  models  are  denoted  by  W(a,c,x,t)  where  a,  c, 
x,  and  t  are  the  Wang  hailstone  parameters  defined  in  the 
previous  chapter.  Parameters  a,  c,  and  t  are  given  in  cm 
and  X  is  dimensionless. 

Scattering  at  3  GHz.  At  3  GHz,  the  wavelength  of  the 
incident  wave  is  10  cm.  The  relative  dielectric  constants 
used  for  ice  and  water  at  3  GHz  are  3.1698  +  JO. 003788  and 
79.6918  +  j25.1975  respectively  (12:1837-1840).  In  this 
section,  scattering  characteristics  of  the  nonmirror 
symmetric  nonlayered  and  layered  Wang  hailstone  models  are 
compared  to  the  scattering  characteristics  of  corresponding 
equivolume  spheres.  Effects  of  different  layer  thicknesses 
and  scatterer  sizes  are  also  presented. 

Figure  14  shows  Wang  model  W(0. 3,0. 6,2,0)  and  an 
equivolume  sphere.  Figure  15  shows  normalized  DSCS  vs. 
scattering  angle  for  the  two  hailstone  models  at  horizontal 
polarization.  Normalization  is  with  respect  to  nc^.  The 
largest  scattering  differences  occur  in  the  forward  and 
backscattered  directions  where  the  DSCS  of  the  Wang  model  is 
about  1.5  dB  below  that  of  the  sphere. 

Figure  16  shows  Wang  model  W(0. 3, 0.6, 2, 0.05).  Figure 
17  shows  normalized  DSCS  vs.  scattering  angle  for 
W(0. 3,0.6, 2,0.05)  and  an  equivolume  layered  sphere  at 
horizontal  polarization.  Uniform  water  layer  thickness  for 
the  equivolume  sphere  is  approximately  0.05  cm.  The  DSCS  is 
greater  for  the  layered  models  than  for  the  homogeneous 
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models  in  Figure  15.  Forward  scattering  differences  are 
more  significant  than  in  the  backscattered  direction.  The 
minimum  return  from  the  Wang  model  is  about  10  dB  greater 
and  10  degrees  to  the  right  of  the  equivolume  sphere's 
minimum  return. 

Figure  18  shows  DSCS  vs.  scattering  angle  at  horizontal 
polarization  for  three  different  hailstone  models.  The 
external  dimensions  of  the  models  are  the  same  while  water 
layer  thickness  varies  from  0  to  0.02  to  0.05  cm.  As  water 
thickness  increases,  DSCS  increases  and  the  minimum  DSCS 
moves  in  the  backscattered  direction.  The  forward 
scattering  characteristics  are  more  affected  by  the  layer 
thickness  than  are  the  backscattered  characteristics. 


SCATTWING  ANGLE  CO60D 


Figure  18.  Normalized  DSCS  vs.  Scattering  Angle  for 
W(0.3,0.6,2,t) 


Figure  19  shows  Wang  models  W( 0 . 1 5 , 0 . 3 , 2 , 0 )  , 


W(0.3,0.6,2,0),  and  W(0.45,0.9,2,0).  The  general  shape  of 
the  three  nonlayered  models  remains  constant  as  their 
physical  dimensions  increase.  Figure  20  shows  DSCS  vs. 
scattering  angle  at  horizontal  polarization  as  scatterer 
size  varies.  DSCS  increases  with  increased  scatterer  size. 

All  plots  in  the  above  section  were  for  a  horizontally 
polarized  incident  wave.  Although  not  shown,  the  vertically 
polarized  scattering  characteristics  were  essentially  the 
same  as  the  horizontally  polarized  scattering 
characteristics . 


Figure  19.  Wang  Models  W ( 0 . 1 5 , 0 . 3 , 2 , 0 ) ,  W ( 0 . 3 , 0 . 6 , 2 , 0 )  and 
W( 0 . 45 , 0 . 9 , 2 , 0 ) 
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Figure  20.  Normalized  DSCS  vs.  Scattering  Angle  for 
W(0. 15,0.3,2,0) ,  W(0.3»0. 6,2,0)  and 
W(0.45,0.9,2,0) 


Scattering  at  30  GHz.  At  30  GHz,  the  wavelength  of  the 
incident  wave  is  1  cm.  The  relative  dielectric  constants 
used  for  ice  and  water  at  30  GHz  are  3.1684  +  jO. 0007352  and 
11.9194  +  j22.5613  respectively  (12:1837-1840).  In  this 
section,  scattering  characteristics  of  the  nonmirror 
symmetric  nonlayered  and  layered  Wang  hailstone  models  are 
compared  to  the  scattering  characteristics  of  corresponding 
equivolume  spheres.  Effects  of  different  layer  thicknesses 
and  scatterer  sizes  are  also  presented.  Forward  and 
backward  scattering  quantities  are  plotted  with  respect  to 
angle  of  incidence.  Normalization  of  all  scattering 
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quantities  are  with  respect  to  nc^. 

Figure  21  shows  Wang  model  W(0.45,0.9,4,0)  and  an 
equivolume  sphere.  Figure  22  plots  normalized  D3CS  vs. 
scattering  angle  for  W( 0 . 45 , 0 . 9 , 4 , 0 )  and  an  equivolume  ice 
sphere  at  horizontal  polarization.  Water  layer  thickness 
for  the  equivolume  sphere  is  approximately  0.05  cm.  DSCS 
varies  greatly  over  the  range  of  scattering  angles. 
Backscattered  DSCS  is  affected  more  by  the  scatterer  shape 
than  the  forward  DSCS. 

Figure  23  shows  Wang  model  W ( 0 . 45 , 0 . 9 , 4 , 0 . 05 ) .  Figure 
24  shows  normalized  DSCS  vs.  scattering  angle  for 
W( 0 . 45 , 0 . 9 , 4 , 0 . 05 )  and  an  equivolume  sphere  at  horizontal 
polarization.  The  DSCS  is  very  different  for  all  scattering 
angles.  Backscattered  DSCS  differs  by  18  dB  and  forward 
DSCS  differs  by  8  dB. 


Figure  21.  Wang  Model  W( 0 . 45 , 0 . 9 , 4 , 0 )  and  Equivolume  Sphere 


SCATTERING  ANGLE  COeflD 


Figure  24.  Normalized  DSCS  vs.  Scattering  Angle  for 
W( 0 . 45 , 0 . 9 , 4 , 0 . 05 )  and  Equivolume  Sphere 


Figure  25  shows  Wang  model  W( 0 . 45 , 0 . 9 , 2 , 0 . 02 ) . 

Figures  26  and  27  show  normalized  DSCS  vs.  scattering  angle 
for  W( 0 . 45 , 0 . 9 , 2 , t )  for  t  =  0,  0.02,  and  0.05  at  horizontal 
and  vertical  polarizations  respectively.  The  forward  DSCS 
does  not  vary  significantly  as  water  thickness  changes. 
Backward  DSCS  varies  greatly  with  layer  thickness.  Forward 
scattering  is  similar  for  horizontal  and  vertical 
polarizations.  Except  for  the  forward  direction,  horizontal 
and  vertical  polarizations  bear  little  resemblance. 


SCATTERING  ANGLE  COeOD 


Figure  27.  Normalized  DSCS  vs.  Scattering  Angle  for 
W( 0 . 45 , 0 . 9 , 2 , t )  Vertical  Polarization 

Figures  28  and  29  show  normalized  DSCS  vs.  angle  of 
incidence  for  W( 0 . 45 , 0 . 9 , 2 , t )  in  the  backward  and  forward 
directions  respectively.  The  incident  wave  impinges  on  the 
top  of  the  hailstone  for  an  angle  of  incidence  equal  to  0 
degrees.  At  90  degrees,  the  incident  wave  approaches  the 
hailstone’s  broadside.  The  incident  wave  impinges  on  the 
bottom  of  the  hailstone  model  at  180  degrees.  Backscattered 
DSCS  is  varies  as  much  as  28  dB  depending  on  the  angle  of 
incidence  and  layer  thickness,  t.  Forward  DSCS  is  much  less 
sensitive  and  varies  only  about  2.5  dB. 
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The  scattering  characteristics  in  the  next  five  graphs 
are  presented  as  functions  of  scatterer  size.  Hailstone 
parameter  c  increases  from  0.3  to  0.9  cm;  parameter  a 
simultaneously  increases  from  0.15  to  0.45  cm.  The  general 
shape  of  the  hailstone  models  do  not  change  as  size 
increases . 

Figures  30  and  31  show  how  forward  and  backward  DSCS 
changes  with  increasing  scatterer  size.  The  forward  DSCS  of 
the  nonlayered  hailstone  model  is  2-3  dB  greater  than  the 
layered  model  over  most  of  the  size  range.  Backward  scatter 
increases  with  increased  scatterer  size  for  the  nonlayered 
hailstone  model.  The  layered  model’s  backward  DSCS  exhibits 
a  lobed  behavior  with  increased  scatterer  size. 
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Figure  30.  Forward  Normalized  DSCS  vs.  Scatterer  Size  for 
W(c/2,c,4,t)  Horizontal  Polarization 


Figure  31.  Backward  Normalized  DSCS  vs.  Scatterer  Size 
for  W(c/2,c,4,t)  Horizontal  Polarization 

Figures  32,  33  and  34  show  scattering,  absorption,  and 
extinction  cross  sections,  respectively,  vs.  scatterer  size 
for  horizontal  polarization.  Extinction  and  scattering 
cross  sections  for  the  layered  model  are  generally  less  than 
for  the  nonlayered  hailstone  model.  Absorption  cross 
section  is  greater  for  the  layered  (ice  and  water )model  than 
for  the  ice  only  nonlayered  model.  Absorption  is  greater 
for  the  layered  model  due  to  the  large  imaginary  part  of 
water’s  relative  dielectric  constant. 
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Figure  34.  Normalized  Extinction  Cross  Section  vs. 

Scatterer  Size  for  W(c/2,c,4,t)  Horizontal 
Polarization 


This  chapter  presented  a  sample  of  the  scattering 
characteristics  that  can  be  calculated  using  the  modified 
EBCM  computer  program.  Scattering  characteristics  vary  with 
incident  wave  frequency;  hailstone  model  size,  shape,  and 
structure;  observation  angle;  and  angle  of  incidence  of  the 
illuminating  wave.  The  scattering  results  are  summarized 
and  recommendations  for  future  work  are  presented  in  the 
final  chapter. 


V.  Conclusions  and  Recommendations 

Introduction 

This  chapter  summarizes  the  results  of  the  preceding 
analysis  and  presents  recommendations  for  future  work  in 
this  area. 

Conclusions 

The  EBCM  was  found  to  be  well  suited  in  determining  the 
scattering  characteristics  of  the  Wang  hailstone  models 
considered.  Conclusions  are  presented  for  each  of  the  two 
incident  wave  frequencies  examined. 

Scattering  at  3  GHz.  Rayleigh  scattering  occurs  at  3 
GHz  where  wavelength  is  much  larger  than  the  hailstone 
models  analyzed.  The  wavelength  is  10  cm  and  the  hailstone 
models  varied  in  size  from  0.6  to  1.8  cm  in  length  along  the 
axis  of  symmetry.  Little  phase  variation  of  the  incident 
wave  occurs  over  the  hailstone  body.  Hailstone  shape  was 
not  a  significant  scattering  parameter.  Horizontal  and 
vertical  polarizations  produced  very  similar  results  for  the 
hailstone  models  examined  at  3  GHz. 

Different  water  layer  thicknesses  of  the  hailstone 
models  produced  significant  forward  and  backward  scattering 
characteristics.  At  3  GHz  the  real  part  of  dielectric 
constant  of  water,  79.6918,  is  much  greater  than  that  for 
ice,  3.1698.  Thus,  the  impedance  mismatch  at  the  hailstone 
surface  is  much  greater  for  the  layered  models  than  for  the 
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ice  only  models.  Forward  and  backward  scattering  is  greater 
for  the  layered  models  than  for  the  ice  only  models  and 
increases  with  increased  water  thickness.  Scattering  at 
other  scattering  angles  was  also  significantly  affected  by 
the  water  layer  and  thickness.  The  bistatic  scattering 
minimum  which  occurred  at  90  degrees  for  the  nonlayered 
hailstone  models  moved  in  the  backscattered  direction  as  a 
water  layer  was  introduced.  The  minimum  moved  further  in 
the  backscattered  direction  as  water  thickness  is  increased. 

The  shape  of  the  nonlayered  hailstone  models  did  not 
affect  scattering  properties  significantly.  The  Wang  model 
and  the  equivolume  ice  sphere  produced  similar  scattered 
fields.  Hailstone  shape  was  more  important  for  the  layered 
models.  Scattering  by  the  layered  Wang  model  was 
significantly  different  than  the  scattering  by  the  layered 
equivolume  sphere.  Scattering  differences  were  most 
pronounced  in  the  forward  scattering  direction. 

Increasing  nonlayered  scatterer  size  produced  larger 
scattered  fields  for  all  scattering  angles. 

Scattering  at  30  GHz.  Resonant  scattering  occurred  at 
30  GHz  where  wavelength  and  scatterer  size  were  comparable. 
The  phase  of  the  incident  radiation  changes  significantly 
over  the  hailstone  length.  The  phase  variations  cause 
complex  interactions  between  the  different  parts  of  the 
hailstones  all  of  which  must  be  considered  in  determining 
the  scattered  fields.  Scattering  characteristics  are  very 
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dependent  on  hailstone  shape,  size  and  structure.  The 
resonant  affects  were  evident  is  the  scattering  calculations 
for  the  different  hailstone  models.  Scattering 
characteristics  at  30  GHz  were  very  dependant  on  incident 
wave  polarization. 

Scatterer  shape  was  more  important  at  30  GHz  than  at  3 
GHz.  A  nonlayered,  nonmirror  symmetric  Wang  model  scattered 
the  incident  wave  much  differently  than  a  corresponding 
equivolume  nonlayered  sphere.  Differences  were  even  more 
pronounced  for  the  layered  Wang  model  as  compared  to  a 
equivolume  layered  sphere. 

Scattering  characteristics  of  Wang  models  changed 
significantly  as  a  water  layer  was  introduced  and  as  water 
thickness  was  changed.  Forward  scattering  was  greater  for 
the  layered  models  than  for  the  nonlayered  models  and 
generally  increased  with  increased  water  thickness.  As  with 
scattering  at  3  GHz,  the  bistatic  scattering  minimums  moved 
in  the  backscattered  direction  as  water  layer  thickness 
increased.  Backward  scattering  changed  significantly  as  the 
angle  of  incidence  of  the  illuminating  wave  changed  from  the 
top  to  the  bottom  of  the  hailstone  models.  Forward 
scattering  changed  only  slightly  with  angle  of  incidence. 
Again  water  thickness  also  affected  scattering  as  angle  of 
incidence  varied. 

The  nonlayered  Wang  models  generally  scattered  more  of 
the  incident  energy  than  the  layered  models  as  scatterer 
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size  increased.  The  layered  models  absorbed  more  energy 
than  the  nonlayered  models  due  to  the  water  layer’s  larger 
complex  dielectric  constant. 

Hailstone  model  size,  shape  and  structure  greatly 
affect  scattering  characteristics  at  30  GHz.  Angle  of 
incidence  is  also  important  particularly  in  backscattered 
directions . 

Recommendations 

Further  research  is  necessary  to  determine  the  actual 
size,  shape  and  composition  of  hailstones.  The  hailstone 
models  examined  in  this  thesis  were  constructed  from  data 
from  one  hailstorm  in  Colorado.  A  larger  sample  size  may 
produce  a  more  accurate  hailstone  model.  Also,  the 
distribution  of  water  on  the  surface  of  the  hailstones  needs 
to  be  examined.  The  water  thickness,  if  present,  needs  to 
be  accurately  modeled  to  obtain  valid  scattering 
calculations . 

The  distribution  of  hailstones  along  a  propagating 
signal's  path  needs  to  be  modeled.  This  thesis  determined 
scattering  characteristics  by  individual  hailstones.  This 
knowledge  along  with  a  hailstone  distribution  model  would 
enable  calculation  of  total  scattering  effects  of  hailstorms 
on  propagating  signals. 

Research  should  be  conducted  to  determine  the 
usefulness  of  the  EBCM  for  predicting  the  radar  cross 
section  of  targets  of  interest  to  the  Air  Force.  The  EBCM 


Modeled  as  Oblate  Spheroids  and  Sphere-Cone-Spheroids 


This  appendix  contains  differential  reflectivity 
calculations  for  a  sample  of  hailstone  models  used  by  Aydin 
and  Seliga.  Aydin  and  Seliga  used  the  EBCM  to  determine  the 
differential  reflectivity  of  hailstones  modeled  as  oblate 
spheroids  and  sphere-cone-spheroids  at  10  cm  wavelengths 
(1:58-66).  A  portion  of  Aydin  and  Seliga* s  work  was 
duplicated  to  test  the  program  used  in  this  thesis. 

Oblate  Spheroid. 

The  differential  reflectivity  (ZDp)  radar  parameter  is 
given  by  Z^jj  =  <*jj/av  where  a^y  are  t^ie  backscatter  cross 
sections  at  horizontal  and  vertical  polarizations, 
respectively.  Z DR  was  calculated  for  various  sizes  of 
layered  and  nonlayered  oblate  spheroids.  The  oblate 
spheroidal  hailstone  model  is  shown  in  Figure  35. 
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Figure  35.  Oblate  Spheroidal  Hailstone  Model 


Figure  36  shows  Zpp  vs.  Equivolume  Diameter  (D)  for  layered 
and  nonlayered  oblate  spheroidal  hailstone  models.  D  is  the 
diameter  of  an  equivolume  sphere  and  is  given  by 
D  =  2a(b/a)^“*  where  a  and  b  are  semiaxes  of  the  spheroid. 
Water  thickness  (t)  is  0  mm  for  the  nonlayered  model  and  0.5 
mm  for  the  layered  model.  The  axial  ratios  of  the  oblate 
spheroids  are  b/a  =  b'/a’  =  0.9.  The  relative  dielectric 
constants  used  for  ice  and  water  are  €r  =  3.1698  +  jO. 003788 
and  79.6918  +  j25.1975,  respectively. 
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Figure  36.  Znp  vs.  Equivolume  Diameter  for  Oblate  Spheroid 


Figure  37  shows  the  sphere-cone-sphere  hailstone  model. 


Parameter  a  is  the  radius  of  the  base  sphere.  Parameter  b 
is  the  radius  of  the  apex  sphere.  Cone  angle  is  given  by  a. 


Figure  37.  Sphere-Cone-Sphere  Hailstone  Model 


ZDR  was  determined  for  various  sizes  of  sphere-cone- 
spheroidal  hailstone  models  and  is  shown  in  Figure  38. 
Equivolume  spherical  diameter  is  given  by  D  =  2.07a  and  is 
varied  from  2  through  26  mm.  The  hailstone  models  for 
Figure  38  have  cone  angle  a  s  60  degrees  and  radius  ratio 
b/a  =  0.5.  The  relative  dielectric  constant  of  the  ice  is 
the  same  as  for  the  oblate  spheroid. 
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Figure  38.  ZDp  vs.  Equivolume  Diameter  for  Sphere-Cone- 
Sphere 

Sphere-Cone-Oblate  Spheroid. 

Figure  39  shows  the  sphere-cone-oblate  spheroid 
hailstone  model.  Parameter  c  is  the  radius  of  the  apex 
sphere.  Parameters  a  and  b  are  the  oblate  spheroid  semiaxis 
lengths.  Cone  angle  is  given  by  angle  a.  Figure  40  shows 
ZDp  vs.  equivolume  spherical  diameter  (D).  The  following 
relationships  hold  for  the  hailstone  models  in  Figure  40: 
c/a  =  0.2,  b/a  =  0.5,  «  =  70  degrees,  and  D  =  1.672a.  The 
relative  dielectric  of  the  model  is  €r 


5.33  +  j  0 . 5  7  . 


The  scattering  calculations  for  the  hailstones  in  this 


appendix  appeared  to  be  exactly  the  same  as  those  published 
by  Aydin  and  Seliga  (1:58-66). 


.  J 
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Appendix  B:  EBCM  FORTRAN  Code  Description 


Introduction 

This  appendix  briefly  describes  the  three  interactive 
FORTRAN  programs  used  in  this  thesis.  Major  function, 
input,  and  output  associated  with  each  of  the  three  programs 
are  presented.  MCON,  MLAYRC ,  and  MAINOR  are  the  three 
programs  discussed. 

Figure  41  shows  the  coordinate  system  used.  The 
incident  wave  propagates  in  the  +z-direction  and  its 
polarization  vector  is  specified  by  amplitude,  A,  and 
angle,  p,  from  the  x-axis  in  the  xy-plane.  Scatterer 
orientation  is  specified  by  angles  ®  and  <t>  where  S  is  the 
axis  of  symmetry  of  the  scatterer.  The  scattered  field  is 
evaluated  in  the  xz-plane  with  0  degrees  in  the  forward  or 


Figure  41.  Coordinate  System 


+z-direction  and  180  degrees  in  the  backscattered  or 
-z-direction.  The  scattering  angle  is  equal  to  90  degrees 
at  the  +x-axis. 

MCON 

MCON  calculates  the  T-matrix  for  a  nonlayered  Wang 
hailstone  model  at  the  incident  wave  frequency  of  interest. 
Gaussian  quadrature  is  used  to  perform  the  required 
numerical  integration.  Bessel  functions  are  computed  using 
downward  recursion.  Legendre  functions  are  computed  using 
upward  recursion.  Matrix  inversions  are  accomplished  using 
a  Gauss- Jordan  inversion  technique. 

Input .  The  interactive  program  prompts  the  user  to 

enter  the  following  variables: 

A:  Wang  hailstone  parameter  (a)  in  cm 

C:  Wang  hailstone  parameter  (c)  in  cm 

LMDA:  Wang  hailstone  parameter  (x) 

WVLTH:  Incident  wavelength  in  cm 

NRANK:  Summation  truncation  value  (N) 

NTHETA:  Number  of  nodes  for  Gaussian  integration 

(multiple  of  2) 

IOUT:  Convergence  testing  code  (0  =  NTHETA,  1  =  NRANK, 

2  =  NM) 

DCNR:  Real  part  of  hailstone  dielectric  constant 

DCNI :  Imaginary  part  of  hailstone  dielectric  constant 

Convergence  Testing.  The  following  procedure  was 
carefully  followed  for  variables  NRANK,  NTHETA,  and  IOUT. 

1.  Set  NRANK  =  10,  NTHETA  =  20,  and  IOUT  =  0. 


The  program  checks  for  convergence  over  NTHETA.  Convergence 
occurs  if  the  DSCS  at  ten  scattering  angles  does  not  change 
by  more  than  one  percent  as  NTHETA  increases  to  its 
specified  value.  If  the  program  responds  with 
***  SOLUTION  HAS  CONVERGED  *** 
go  to  Step  2.  If  convergence  is  not  achieved  run  the 
program  again  with  NRANK  =  10,  IOUT  =  0,  and  NTHETA 
increased  by  2. 

2.  Set  NRANK  =  10,  IOUT  =  1  and  NTHETA  equal  to 
its  convergent  value  from  Step  1.  The  program  will  run  and 
check  for  convergence  over  NRANK.  Convergence  occurs  if  the 
DSCS  at  ten  scattering  angles  does  not  change  by  more  than 
one  percent  as  NRANK  increases  to  its  specified  value.  If 
the  program  responds  with 

***  SOLUTION  HAS  CONVERGED  *** 

go  to  Step  3.  If  convergence  is  not  achieved  increase  NRANK 
and  run  the  program  again  with  IOUT  =  1  and  NTHETA  equal  to 
its  convergent  value  from  Step  1. 

3.  Set  NRANK  and  NTHETA  equal  to  their  convergent 
values  and  IOUT  =  2.  The  program  will  run  to  convergence 
over  NM.  Variable  NM  corresponds  to  index  m  of  the 
associated  Legendre  functions. 

Output .  The  output  of  MCON  consists  of  two  data  files: 
EBCM.DAT  and  T.DAT.  EBCM.DAT  contains  the  input  variables 
used  in  computing  the  T-matrix.  T.DAT  contains  the  computed 


T-matrix . 


MLAYRC 


MLAYRC  computes  the  T-matrix  for  a  layered  Wang 
hailstone  model  at  the  incident  wave  frequency  of  interest. 
Numerical  techniques  used  in  MLAYRC  are  the  same  as  those 
used  in  MCON. 

Input .  The  interactive  program  prompts  the  user  to 
enter  the  same  variables  as  in  program  MCON.  Wang  hailstone 
parameters  a,c,  and  x  correspond  to  the  external  dimensions 
of  the  layered  hailstone  model.  MLAYRC  prompts  the  user  to 
enter  water  layer  thickness  (t)  and  the  real  and  imaginary 
parts  of  the  dielectric  constants  for  each  layer  of  the 
hailstone  model. 

Convergence  Testing.  Convergence  testing  for  MLAYRC  is 
accomplished  using  the  same  procedure  as  in  MCON. 

Output .  Same  as  for  MCON. 

MAI NOR 

MAINOR  computes  DSCS ,  scattering,  absorption,  and 
extinction  cross  sections  for  a  specified  incident  wave  and 
hailstone  orientation.  All  scattering  quantities  are 
normalized  to  where  c  is  the  Wang  hailstone  parameter. 

Input .  MAINOR  uses  the  output  files  from  MCON  and 
MLAYRC  to  calculate  the  required  scattering  quantities. 
MAINOR  also  interactively  prompts  the  user  to  enter  the 
following  variables: 

AMPLAB :  Amplitude  (A  from  Figure  41)  of  incident  wave 

ANGLAB:  Polarization  angle  ( v>  from  Figure  41)  of 

incident  wave 


‘j?£j 


M  '  .  "I1 


IOUT:  0  (no  output  file)  or  1  (output  file  SCAT.DAT) 

DLTANG:  Scattering  angle  increment 

NUANG:  Number  of  scattering  angles 

NCASES:  number  of  hailstone  orientations 

ANGINT:  Hailstone  orientation  angle  (©  from  Figure  41) 
ANGINP:  Hailstone  orientation  angle  («  from  Figure  41) 
MAINOR  will  prompt  the  user  for  ANGINT  and  ANGINP  for  each 
desired  hailstone  orientation. 

Output .  Normalized  DSCS  (for  each  specified  scattering 
angle),  scattering,  absorption,  and  extinction  cross 
sections  are  displayed  on  the  user’s  terminal  screen.  If 
requested,  DSCS  calculations  are  written  to  data  file 
SCAT.DAT. 
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